Reducing cadmium bioaccumulation in Theobroma cacao using biochar: basis for scaling-up to field

The intake of Cd-enriched food is the main Cd pathway for the nonsmoking population. In some cases, Cd bioaccumulates in edible plant parts which comprise risk to consumers, because of Cd is a harmful heavy metal that can cause potent environmental and health hazards. For instance, Cd enrichment of cacao seeds have led to Cd enrichment of cacao-based products. In Latin America and the Caribbean, Cd bioaccumulation in cacao seeds occurs in different regions with diverse edaphoclimatic conditions, which makes it difficult to select soil remediation alternatives. Limited resources require that potential amendments must be carefully investigated through laboratory and/or greenhouse conditions before scaling up to field experiments. In this study, we evaluated the effectiveness of four biochars: coffee-, quinoa-, and inoculated- and palm-biochar, derived from three feedstocks: coffee husk, quinoa straw, and oil palm residues, respectively. Biochars were applied in two rates (1 and 2% w/w) in two soils, one moderately acidic and one slightly alkaline (Cd-spiked and non-spiked). CCN-51 cacao plants were used for the greenhouse experiment. After 130 days, biometric parameters, the bioavailability of Cd in the soil, and the concentration of Cd and mineral nutrients in the plants were measured. Quinoa biochar at the 2% significantly decreased (P < 0.01), by ∼71%, bioavailable Cd in moderately acidic and slightly alkaline soils, and leaf-Cd by ∼48%. Soil pH, electrical conductivity, and effective cation exchange capacity were significantly (P < 0.01) correlated with bioavailable soil and leaf-Cd. Biochar characteristics, such as ash contents, basic cations content, and surface functional groups could be used as indicators for the selection of biochars to reduce Cd uptake by cacao. Additionally, application of quinoa derived biochar provided P and K, which could increase productivity to offset mitigation costs. Overall, incorporation of quinoa biochar at 2% rate is effective for lowering bioavailable Cd in different soil types which reduces leaf-Cd in cacao plants.


Introduction
Cadmium (Cd) is a natural-occurring soil metal that can potentially affect human health (Smolders and Mertens, 2013). Chronic exposures of Cd can cause bone disease, endocrine disruption and has been shown to be a potential carcinogen (Pan et al., 2009;Suhani et al., 2021). Although there are different Cd exposure routes (such as inhalation and cigarette smoking), the intake of Cd-enriched food is the main Cd pathway for the nonsmoking population (WHO, 2011;European Food Safety Authority, 2012;Suhani et al., 2021). This phenomena have been reported in edible plant parts of several crops like: wheat grains (Rizwan et al., 2016), rice grains , and cacao seeds (Argüello et al., 2019), the mainstay for the chocolate industry. Environmental authorities have established concentration limits for Cd in chocolate, cocoa powder, and derivatives. For instance, the European Union and the State of California (USA) have established limits between 0.1 -0.8 and 0.4-0.96 mg kg À1 Cd in final product, respectively (European Commission, 2014;Maddela et al., 2020). These limits for Cd are related to the final product, not to the concentration in the cacao seeds, however traders commonly applied levels of Cd in the beans (Bravo et al., 2022). For cacao seeds, commonly applied threshold for export to the EU market is 0.6 mg kg À1 , an unofficial industry threshold (Argüello et al., 2019;Vanderschueren et al., 2022). In Latin America and the Caribbean the Cd concentrations in cacao seeds averaged 0.9, 1.0, 2.5, and 4.0 mg kg À1 in Ecuador (Argüello et al., 2019), Peru (Ar evalo- Gardini et al., 2017), Honduras (Gramlich et al., 2018), and Colombia (Bravo et al., 2022), respectively. These values are much higher than African countries in which Cd concentration averaged 0.3 mg kg À1 (Bertoldi et al., 2016). In this context, Cd concentration in cacao seeds is a menace for the sustainability of small-scale farmers (Meter et al., 2019). Cadmium in cacao seeds is strongly related to Cd concentration in soils, key soil properties (e.g. soil pH, soil organic matter, and P content) (Chavez et al., 2016b;Argüello et al., 2019;Bravo et al., 2022), and content of other elements like Zn (Argüello et al., 2019). Soil Zn content may favor or decrease Cd bioaccumulation depending on Zn bioavailability and Zn/Cd ratio, since the Cd transporters are nor specific, and the uptake of this element has been related to Zn transporters (Argüello et al., 2020).
Then, in-situ Cd immobilization in cacao-growing soils using amendments, has been proposed as a promising option (Chavez et al., 2016a;Argüello et al., 2019); nonetheless, this alternative has been poorly explored to date. A widely known amendment to control soil contamination is biochar (Brewer et al., 2011). The immobilization mechanisms and their effectiveness depend on several factors e.g. soil pH (Xiao et al., 2019) and material characteristics (ash, basic cation, and P content) (L opez et al., 2020;Houssou et al., 2022). Soil pH is among the most critical parameters influencing the effectiveness of biochar in Cd mitigation (Xiao et al., 2019). Uchimiya et al. (2011) found that cation exchange mechanism played a dominant role when biochar was used to remediate acidic soils. Likewise, the contribution of alkaline substances by biochars plays a vital role when biochar is applied to alkaline soils. However, less information was available for the method of improving Cd immobilization efficiency taking into account the pH of the targeted soil (Jia et al., 2022). Cacao production areas cover a wide diversity of soils with contrasting properties (e.g. acidic and neutral to alkaline soils) (Chavez et al., 2016a;Argüello et al., 2019), which represents a challenge for the selection of soil amendments to reduce plant-available soil Cd. There is little information about Cd mitigation in cacao-growing soils, especially in neutral to alkaline soils. Only one previous study evaluated the effect of biochar application on Cd uptake by cacao plants (Ramtahal et al., 2019).
In order to better understand the complex interplay between Cd, soil and biochar, in a broader scope, there is a need to conduct systematic studies evaluating different types of biochars at various rates under distinctive soil conditions. Developing criteria for selecting amendments for field studies based on the chemical properties of the material and/or soil is urgent. The aim of this work was to evaluate the effects of four biochars at two low rates (1 and 2% w/w) on the properties of two distinctive cacao-growing soils (slightly alkaline and moderately acidic) and their relationship with Cd bioavailability and uptake by cacao plants. Such research could help to alleviate excessive Cd bioaccumulation in cacao growing regions.

Soil selection and sampling
For an easier follow-up and understanding of the methodology, an experimental overview can be found in Supplementary information ( Figure S1). Soils were collected at two cacao farms where elevated bean Cd were reported (Argüello et al., 2019). Slightly alkaline soil (pH 7.30) was collected in Manabí province (Northeast) whereas moderately acidic soil (pH 5.52) was collected in Azuay province (Southeast). At these two cacao farms, Cd concentration in cacao bean exceeded the Cd threshold of 0.60 mg kg À1 (the European Union limit for Cd in chocolate), considered excessive for trading. A total of 500 kg of surface soil (0-15 cm) was collected at each farm by composite sampling (twelve subsamples within each farm). Soils were then air dried, disaggregated, sieved (2-mm mesh) and mixed to obtain a homogenous substrate. Soil pH and electrical conductivity (EC) were measured using 2:1 water:soil ratio. Soil organic matter (SOM) was determined using an elemental analyzer (Elementar Analyzer, Germany). Pseudo-total concentration of Cd and mineral nutrients were measured by Inductively Coupled Plasma Optical Emission Spectrometry (ICP-OES, Optima 5300 DV, PerkinElmer, USA) after digesting 0.25 g of soil sample with 10 mL of Aqua Regia (1:3 HCl:HNO 3 ) at 90 C for 1 h and 140 C for 3 h.

Biochars
The four biochars used in this study were produced by Bioenergía de los Andes (Ecuador) using a prototype of a modular auger reactor for the pyrolysis process, integrated with a combustion reactor . Biochars were produced from three feedstocks: coffee husk (Coffee-BC, 550 C), oil palm residue (Inoculated-BC, 600 C and Palm-BC, 600 C), and quinoa straw (Quinoa-BC, 550 C). Inoculated-BC is a commercial mixture of oil palm residue biochar with a microbiology consortium. More details about the production of these biochars were reported in L opez et al. (2020). Biochar pH and EC were determined using 1:10 soil to water ratio after 60 min. Ash content was determined by loss on ignition. Carbon (C) and nitrogen (N) content were determined using an elemental analyzer. Pseudo-total concentration of Cd and mineral nutrients in extracts were measured by using ICP-OES as indicated above. Fourier-transform infrared spectroscopy analysis (FTIR) was used to characterize the biochar active sites, for which an infrared spectrometer (PerkinElmer, model spectrum Two V10.4.2) equipped with an Attenuated Total Reflection (ATR) accessory (PerkinElmer) was employed, operating in the spectral range 400-4000 cm À1 with a resolution of 4 cm À1 .

Soil preparation and incorporation of biochars
Before biochar incorporation, the moisture contents of soils were kept at 40% water holding capacity (WHC) with deionized water (DW). Then, each biochar was incorporated in the soils at two rates (1 and 2% dry w/ w) and soils without biochar were kept as controls. Later, soil moisture was adjusted to 80% of WHC with DW and equilibrated for 7 days. After equilibrium, the soils with and without biochars were divided in two parts, one part was kept as it is in the field (natural Cd concentration) and the other part was spiked with Cd (enriched conditions) using 1.00 mg Cd [from Cd(NO 3 ) 2 4H 2 O] per kg of soil. The former to evaluate the effectiveness of biochars under "natural" conditions and the latter to determine the potential of biochars to mitigate "enriched" conditions. The Cd concentration in spike solution was selected based on total soil Cd reported in Argüello et al. (2019) and was calculated to include the 97.5 percentile of the soil Cd reported in Ecuador, approximately 1.40 mg kg À1 . The spiking solution was spread onto the soils while mixing thoroughly using a plastic spatula. The pseudototal Cd concentration in slightly alkaline soil and moderately acidic soil were 0.57 (AE0.06) and 0.61 (AE0.12) mg kg À1 , respectively. After spiking, average pseudototal Cd concentration in slightly alkaline soil and moderately acidic soil samples were 1.38 (AE0.44) and 1.64 (AE0.14) mg kg À1 , respectively, reaching the targeted concentrations. All treatments were brought to 80% WHC with DW and allowed to further equilibrate for 8 days.

Greenhouse experiment
After the 8-day equilibrium, 4 kg of soil of each treatment was transferred into polyethylene bags (30 Â 40 cm width and depth, respectively) and then 6-month cacao plants were transplanted into each bag. Plant material (cv. CCN-51) was propagated by rooted cutting to avoid the intrinsic genetic variability of seed propagation in cacao. 108 experimental units were arranged in a completely randomized design with three replicates. During the greenhouse experiments, plants were watered (with DW) periodically to maintain 80% WHC. Cacao plants were fertilized monthly using 1.76 g bag À1 of Multi K (13% N, 0% P, 14% K) and 0.44 g bag À1 of Rafos (12% N, 24% P, 12% K, 2% Mg). Plants were allowed to grow under greenhouse conditions, 12 h light and 12 h dark, and average day temperature of 30 AE 1 C, for 130 days. During the greenhouse experiments soil subsamples were withdrawn from each bag at 15, 30, 60 and 90 days to determine the pH dynamics.

Plant harvest and analyses
After 130 days, stem length, root length and stem diameter, and shoot and root dry weight were measured. Cacao plants were separated into leaves, stem and roots and subsequently washed. Leaves and stem were washed with DW, whereas roots were washed with tap water and then with DW, and subsequently immersed in a 0.01 M EDTA solution for 5 min followed by rinsing with DW. Roots, stems and leaves were dried at 65 C for 48 h and dry biomass was measured. Dried tissues were ground to determine elemental composition of each plant part. In brief, 300 mg of pulverized dry biomass were digested with 5.0 mL of HNO 3 at 80 C for 6 h. After digestion, samples were diluted with DW and passed through a 0.45 μm filter paper. Mineral nutrients and Cd in the filtrates were determined by ICP-OES.

Post-harvest soil samples and chemical analyses
After harvest, soil samples were collected and pH, EC, SOM, effective cation exchange capacity (CECe), and Cd pools were measured. Soil pH, EC, and SOM were also measured as indicated above for biochars. Soil CECe was determined as a sum of exchangeable base cations extracted with 1 M NH 4 OAc solution. Three different extraction procedures were used to study the Cd pools in soil samples. The extraction methods were selected based on their extraction power and targeted pools, which are explained as follows: 0.01 M CaCl 2 extract water-soluble and readily exchangeable metals (CaCl 2 -Cd), 1 M NH 4 OAc extract solubleexchangeable metals (NH 4 OAc-Cd), and 0.1 M HCl extract exchangeable and acid-soluble metals (HCl-Cd) (Adamo and Zampella, 2008). For the CaCl 2 -, NH 4 OAc-, and HCl-Cd extractions were using 1:5, 1:5, and 1:20 sample to extractant ratio and equilibration time of 2, 2, and 3 h, respectively. Cadmium in supernatants was determined by ICP-OES.

Quality assurance, quality control, and data analysis
To assure the reliability of the results, quality assurance and quality control protocols were included throughout the analytical processes (see Supplementary Description S1). Two-way ANOVA was used to compute interactions between biochar type and rate. To determine the main effect of biochar or rate, one-way ANOVA was employed. Post hoc tests (Tukey and Dunnet) were applied to calculate differences among groups. For all the analyses, statistical differences were set at 5 %. Multivariate analysis between soil properties, soil Cd pools, and leaf-Cd was used to understand the interactions between soils and biochars. All data was analyzed using JMP ® Pro version 13.1.0 software. When necessary, the continuous variables were transformed to comply with statistical assumptions. Reduction factors (Rf) were calculated using the ratio between Cd concentration in control divided by the corresponding concentration in treated soil/plant. Rf aims to illustrate the proportional change caused by the treatment as compared to the control, i.e., values over 1 imply a reduction in Cd concentration whereas values lower than 1 the opposite (Vanderschueren et al., 2021).

Soil and biochars properties
The physicochemical properties of the soils and biochars are listed in Table 1. The two soils presented contrasting characteristics, mainly pH. Biochar pH was similar among all the studied biochars. In contrast, properties such as EC, ash, basic cations, and carbon content, displayed greater variability. Carbon content was highest in Inoculated-BC (67.8 %) and lowest in Quinoa-BC (56.8 %) whereas biochar pH and mineral constituents (e.g. ash and basic elements) were highest in Quinoa-BC and Coffee-BC. The C/N ratio in Quinoa-and Coffee-BC was on average 40% lower compared to Palm-and Inoculated-BC. Table S1 shows the result of the FTIR analysis. A greater number of functional groups were detected on the surface of Quinoa-BC in comparison with Coffee-, Palm-, and Inoculated-BC. Pseudototal Cd concentration was below the limit of detection of ICP-OES in all the studied biochars.

Soil properties and Cd pools
The effect of biochar application and rate on soil pH is shown in Figure 1 ab, Figure S2 and Table S2. All biochars increased the soil pH, however, the liming effect was not maintained for a long period, i.e., >14 days ( Figure S2). Only Quinoa-BC at 2% showed a consistent and more extended liming effect compared to Coffee-, Inoculated-, Palm-BC, and the control, regardless of soil type ( Figure S2). Quinoa-BC at 2% significantly increased (P < 0.05) soil pH after 130 days by~1.11 units and 0.90 units in slightly alkaline and moderately acidic soils, respectively. Quinoa-BC at 2% raised soil pH by 1.51 in non-spiked slightly alkaline soil (Figure 1a), 0.85 in non-spiked moderately acidic and 0.93 in Cdspiked moderately acidic soils (Figure 1b). No statistical difference was detected in soil pH in Cd-spiked slightly alkaline soils ( Figure 1b).
The effect of biochar application and rate on the CECe is shown in Figure 1(c and d) and Table S2. In line with the pH trend, CECe was higher in biochar-amended soils. Quinoa-BC at the 2% rate increased CECe by 13 % in non-spiked slightly alkaline soil, 11 % in non-spiked moderately acidic soil (Figure 1d), and 16 % in Cd-spiked moderately acidic soil (Figure 1d). There was no statistical effect (P > 0.05) in the application of biochars in Cd-spiked slightly alkaline soils ( Figure 1c). Soil EC increased proportionally with application rate (Figure 1 [e and f] and Table S2). Application of Coffee-and Quinoa-BC showed the most significant (P < 0.05) increases in soil EC in all soils with~27.7% higher EC in slightly alkaline soil with (Figure 1e) and~26.1% greater EC in moderately acidic soils ( Figure 1f). Cd-spiked slightly alkaline soil showed a significant increment in EC (P < 0.05) for all biochar-amended soils ( Figure 1e). In Cd-spiked moderately acidic soil this phenomenon was not observed.
The application of biochars significantly increased (P < 0.05) SOM regardless of soil type or Cd treatment ( Figure 1gh and Table S2). On average, SOM increase was in line with the application rate, which is due to the organic carbon inputs from biochar (Table 1). Our results showed SOM increased between 18 and 77% in slightly alkaline soil ( Figure 1g) and between 31 and 82% in moderately acidic soil ( Figure 1h).
Cadmium pools in soils were determined at 130 days using three single-step extraction methods as shown in Table 2. Water-soluble and easily exchangeable Cd (CaCl 2 -Cd) and exchangeable Cd (NH 4 OAc-Cd) represent the directly bioavailable and reactive pools of Cd in the soil, respectively, while HCl-Cd represents the least available pool (Pan et al., 2016). Quinoa-BC at 2% significantly decreased soluble and exchangeable Cd pools in soils by~71% ( Table 2). Effects of biochar application and rate on reduction factor (Rf) of soil Cd pools are shown in Table 3 and  Table S3. Increasing the biochar rate significantly reduced (P < 0.05) CaCl 2 -and NH 4 OAc-Cd. The largest Rf in CaCl 2 -Cd was observed with application of Quinoa-BC at 2% for all soils, under basal and Cd-enriched conditions. The Rf in soils treated with Quinoa-BC at 2% was~4-fold higher than the other biochars. CaCl 2 -Cd concentration in non-spiked slightly alkaline soils were lower than the detection limit of the ICP-OES and therefore are not reported. Similarly, NH 4 OAc-Cd was significantly (P < 0.05) reduced by the application of Quinoa-BC at 2%. The corresponding Rf for this pool was~2-fold higher for all slightly alkaline and moderately acidic soils. HCl-Cd was reduced only in Cd-spiked moderately acidic soil with the application of Quinoa-BC.

Leaf-Cd content
Cadmium concentration in cacao leaves and the corresponding Rf for plant tissue are shown in Table 2 and Figure 2, respectively. Cadmium concentration in stems and roots are presented in Table S4. In general, cacao plants tend to compartmentalize Cd into the aerial parts. In line with the bioavailable Cd pools in soils, Cd concentration in cacao leaves was~1.4-fold higher in non-spiked moderately acidic soil compared with non-spiked slightly alkaline soil, whereas in Cd-spiked moderately acidic soil the value was~2-fold higher compared with Cd-spiked slightly alkaline soil. Overall, CaCl 2 -Cd and NH 4 AOc-Cd best represented the bioavailable fractions since they were best correlated with leaf-Cd ( Figure 3). Effects of biochar application and rate on leaf-Cd concentration and Rf are presented in Table S5, Table S6, Table 2, and Figure 2. The application of Quinoa-BC at 2% significantly (P < 0.05) lowered Cd concentration in leaves in moderately acidic soil, Cd-spiked moderately acidic soil, and Cd-spiked slightly alkaline soil. In contrast, Cd in leaves in non-spiked slightly alkaline soils did not differ statistically but a lower concentration in the Quinoa-BC at 2% treatment was observed (Table 2).
To clearly show the impact of biochar application on leaf-Cd, the Rf was calculated (Figure 2 and Table S6). In line with leaf-Cd concentration, Quinoa-BC at 2% rate increased (P < 0.05) Rf in all treatments. Average Rf values in Quinoa-BC treatment (at the 2% rate) were 2.01 and 1.18 in non-spiked and Cd-spiked slightly alkaline soils, respectively. The corresponding Rf values for non-spiked and Cd-spiked moderately acidic soils were 1.80 and 2.20, respectively. In addition, with the Cd-spiked moderately acidic soil, Palm-BC at 2% showed similar Rf as Quinoa-BC.

Plant growth and mineral nutrition
No significant effects were found for the biometric parameters measured in the plants under the greenhouse conditions evaluated (P > 0.05, Table S8). The application of Coffee-, Quinoa-, and Palm-BC increased (P < 0.05) plant-P by at least 44% and plant-K over 31% in the moderately acidic soil (Table S8) (P and K deficient soil, Table 1). Similarly, Quinoa-BC at 2% increased plant-K by over 20 % in the slightly alkaline soil.

Biochar properties varied according to feedstock materials
The variability found in the properties of biochars is explained by the different feedstocks and production conditions (L opez et al., 2020). The particular characteristics of each biochar suggest that their effect as a soil amendment could be variable, for instance can impact the ability of these materials for Cd immobilization (L opez et al., 2020). Higher C/N values of Palm-and Inoculated-BC indicate greater aromatic compounds compared to Quinoa-and Coffee-BC. This suggests that the latter two may have a higher Cd adsorption capacity. Low aromaticity provides more sorption sites due to cation (Cd) interaction with existing π-electronic systems from C¼C bounds of the aromatic structure of the biochar (L opez et al., 2020). The oxygen-containing functional groups (hydroxyl, carboxylate, and phenolic hydroxyl) detected in Quinoa-BC could play an important role in the Cd adsorption via Cd-ligands coordinations bonds (L opez et al., 2020). On the other hand, a higher ash content is related to a higher alkalinity of the materials (Yuan et al., 2011), which could affect soil pH in the long term such as a significant increase in pH for acidic soils. Moreover, a higher content of basic cations indicates that Quinoa-BC could exchange these alkali elements with Cd in solution thus decreasing bioavailable Cd (Gao et al., 2019;L opez et al., 2020).
Due to the complex nature of Cd chemistry in soils, identifying the mechanisms that affectively bind Cd to biochar to prevent its uptake by cacao plant is an important task. In a previous study we found that one of the main mechanisms for competitive binding from the soil to the plant are mostly due to ion exchange with surface cations and precipitation as insoluble matter (L opez et al., 2020), which are consistent with the strong dependencies of the Cd absorption with the soil pH (Smolders and Mertens, 2013). The results revealed that chemical properties of the soils were significantly affected by the addition of biochars, especially Quinoa-BC. The higher alkali nature, ash and basic cations content, and superficial functional groups of Quinoa-BC (Table 1 and Table S1) may explain its superior behavior to raise and maintain the soil pH comparing to the other biochars. Previous reports have related the liming effect of biochars with the presence of basic cations (Ca 2þ , Mg 2þ , K þ , and Na þ ), carbonates (CO 2-3 and HCO -3 ) and soluble organic compounds (RCOO-and RO-) (Fidel et al., 2017). In our study, we observed a higher increase in soil pH with the incorporation of Quinoa-BC at 2%. With similar application rates, researchers reported increments of 0.05 and 0.25 units in alkaline soils (Azhar et al., 2019; Xiao et al., 2019) and 0.24 and 1.08 units in Figure 1. Changes on soil pH, effective cation exchange capacity (CECe), electrical conductivity (EC), and soil organic matter (SOM) in coffee-, inoculated-, quinoa-, and palm-biochar amended soils at 1 and 2% rate at 130 days after greenhouse experiment. Slightly alkaline soil (a, c, e and, g). Moderately acidic soil (b, d, f and, h). Error bars are standard errors (SE) of the means (n ¼ 3). Asterisks above the bar indicate a significant (P < 0.05) difference between the treatment and the control by Dunnett's test.
acidic soils (Houben et al., 2013;Shen et al., 2016). This confirms a strong and prolonged liming power of Quinoa-BC. For the Cd-spiked alkaline soils, a stable pH behavior (Figure 1b) could be due to the buffering of the alkaline materials, which is supported by the significant (P < 0.05) increase in EC for the Cd-spike alkaline soil (Figure 1e) (Gao et al., 2019).
Application of biochars increased SOM due to the organic carbon inputs from biochar (Shi et al., 2020). We found a 10-fold increase in SOM in moderately acidic soil as compared to the slightly alkaline soil. The application of biochar into moderately acidic soil could ameliorate soil acidity and increase SOM by reducing CO 2 release from soils (Bashir et al., 2018). Similarly, increases in SOM between 17 and 50% have been reported with the application of rice straw, rice hull, and maize stover biochars at application rates between 1.5 and 3% (Bashir et al., 2018).
The large increase in CECe by application of Quinoa-BC is explained by the raised in soil pH and SOM (Figure 1ab and 1gh, respectively), which enhances the negative soil surface charges and the surface functional groups present in the Quinoa-BC (Table S1) (Shi et al., 2017). Moreover, the concentration of basic cations in Quinoa-BC is higher compared to Coffee-, Inoculated-, and Palm-BC (Table 1), suggesting a greater cation exchange capacity (L opez et al., 2020). A higher EC was expected because these two biochars showed a higher content of basic cations as compared to Inoculated-and Palm-BC (Table 1). A release of weakly bound basic cations from the biochar matrix to the soil solution increases the EC values (Khan et al., 2017). The exchange between basic cations and labile Cd 2þ might explain such increase (Gao et al., 2019). In Cd-spiked moderately acidic soil this phenomenon was not observed probably because of the low base cations of this soil (Table 1).

Soil cadmium pools affected by the incorporation of biochars
Changes in soil chemical properties impacted Cd pools and therefore Cd bioavailability in both soils. Soil pH, CECe, and SOM are key properties governing available Cd in soils (Smolders and Mertens, 2013). Our results showed that, after biochar addition, Cd pools, i.e., 0.01 M CaCl 2and 1 M NH 4 OAc-extractable Cd, decreased due to the increase of key soil properties, i.e., soil pH (Table 4). We found a significant (P < 0.01) and a negative relationship between bioavailable pools with soil pH, CECe, exchangeable basic cations, and EC (Table 4). Higher soil pH increases the sorption capacity of Cd in amended soils. This is caused by the Table 2. Leaf-Cd concentration and CaCl 2 -, NH 4 OAc-, and HCl-Cd concentration in non-spiked and Cd-spiked slightly alkaline and moderately acidic soils amended with coffee-, inoculated-, quinoa-, and palm-biochar at 1 and 2% rate at 130 days after greenhouse experiment. Values are averages of three replicates (n ¼ 3) and AESE represent the standard errors of the mean. Different letters indicate statistical significance at P < 0.05 (Tukey's test) of the soil and Cd treatment group. Data of CaCl 2 -Cd for non-spiked slightly alkaline soil below the detection limit of equipment. variable surface charge, especially in moderately acidic soils (Wang and Siu, 2006). According to our results, soil pH is a good predictor of bioavailable Cd (Figure 3). Based on Figure 3, a pH value equal to or higher than 6.0 significantly reduce (up to 80%) the concentration of bioavailable Cd. Thus, biochars capable of increasing soil pH above this value may be good candidates for remediation in moderately acidic soils. A strong correlation between CaCl 2 -extractable Cd, CECe and basic cations (Table 4) supports our hypothesis that ion exchange was a major mechanism in the immobilization of Cd in soils amended with Quinoa-BC. Interestingly, there was significant (P < 0.01) positive correlation between SOM and both CaCl 2 -and NH 4 OAc-Cd (Table 4). This could be explained by the potential release of low weight organic molecules from biochars which could favor the formation of soluble organometallic complexes and increase soil-plant transfer (Hooda, 2010). Previous studies have concluded that biochars with different feedstock and/or pyrolysis conditions are needed in order to reduce bioavailable Cd in acidic and alkaline soils. Qian et al. (2019) found that the immobilization of Cd and Zn in acidic and alkaline soils require the use of different biochars. Authors found that wheat biochar produced at low temperatures (350 C) was more suitable for acidic soils, whereas biochar produced at high temperatures (650 C) were more effective for alkaline soils. Biochar pH and phenolic functional groups were greater in low-temperature biochar which favor Cd immobilization in acidic soil. On the contrary, our results indicated that a biochar produced from quinoa straw, was effective in lowering available Cd in both soils. Uchimiya et al. (2011) found that different properties, such as cation exchange and alkaline substances, are key features in a biochar for its effective use in lowering bioavailable soil Cd.
Our findings suggest that the changes in soluble and exchangeable Cd in Quinoa-BC amended soils are related to direct and indirect immobilization mechanisms. Cadmium sorption could occur directly on the biochar particles by formation of complexes with functional groups present in the amendment (e.g. C-COOCd þ and C-OCd þ ) (Gong et al., 2021). In the other hand, as soil pH and CECe increases, Cd may be sorbed to variable-charge solids (Fe and Mn oxides) under formation of mononuclear, monodentate binding (e.g. FeOH þ Cd 2þ -> FeOHCd þ þ H þ ) and by binuclear, bidentate binding (e.g. 2FeOH þ Cd 2þ -> (FeO) 2 Cd þ 2H þ ) (Borggaard et al., 2019). Additionally, the increase in SOM favors the formation of Cd-ligand complexation with surface functional groups like carboxylic, hydroxylic, and phenolic (e.g. Cd 2þ þ R-OH -> Cd-RO þ H þ ) (Hamid et al., 2020). In the Cd-spiked moderately acidic soil, a significant change in HCl-extractable Cd may imply the formation of Cd-bearing minerals (Bashir et al., 2018). It was observed that high Cd concentration in solution favors the formation of Cd-bearing minerals, such as CdCO 3 . This occurs at the surface of the Quinoa-BC due to the reaction between Cd and carbonatious materials (L opez et al., 2020). Based on these results, properties such as a high base cation and ash content are desirable in biochars for the immobilization of Cd in Table 3. CaCl 2 -, NH 4 OAc-, HCl-Cd reduction factor (Rf) values in non-spiked and Cd-spiked slightly alkaline and moderately acidic soils amended with coffee-, inoculated-, quinoa-, and palm-biochar at 1 and 2% rate at 130 days after greenhouse experiment. Values are averages of three replicates (n ¼ 3) and AESE represent the standard errors of the mean. Different letters indicate statistical significance at P < 0.05 (Tukey's test). Data of CaCl 2 -Cd for non-spiked slightly alkaline soil below the detection limit of equipment.  cacao-growing soils. In addition, the capability of the biochar to maintain a high soil pH for a prolonged period of time is crucial. The latter is associated with the functional groups present in the material, which help to modify the soil buffering capacity (Shi et al., 2017). After the sorption of Cd in the biochar, this potentially toxic element continues in the soil, however, it will be in a less mobile form, reducing its bioavailability. Because immobilization in these materials occurs by mechanisms such as co-precipitation (L opez et al., 2020), Cd is less likely to return to the soil solution, but as mentioned above, this will be conditioned, for instance, by changes in soil pH. It is necessary to evaluate over time, how the Cd sorbed in the biochar could be exchanged in the soil solution, so reapplication of the biochar over time would be required.

Quinoa biochar reduced plant-Cd
A lower bioavailable soil Cd reduced the Cd concentration in cacao plants, specifically in the soil treated with Quinoa-BC at a 2% rate. Since cacao plant binds bioavailable Cd in the soil through the root, where TcNRAMP5 has been identified as an important transporter for Cd uptake from the soil solution (Ullah et al., 2018). Overall, CaCl 2 -Cd and NH 4 AOc-Cd best represented the bioavailable fractions since they were best correlated with leaf-Cd (Figure 4). These extraction procedures are simple to use and commonly apply to determine readily bioavailable Cd in soils (Adamo and Zampella, 2008).
Cadmium content in leaves was highly correlated with Cd concentration in roots and stems ( Figure S3, r > 0.62, P < 0.01), which is explained by the fact that the cacao plant promotes some tolerance strategies such as uptake and transfer of Cd from the root to the shoot (Oliveira et al., 2022). It can be established that cacao is very efficient in translocating Cd to the leaves. Thus, avoiding the uptake of Cd by the root will decrease the concentration in stems and leaves. In addition, it has been reported that the differential Cd content in the cacao plant can be attributed to the preferential binding of Cd to compounds such as phytate, organic acids, and polyphenols (Vanderschueren et al., 2022).
The discussion will focus on leaves rather than stem or roots. Our results indicate that at a rate of 2%, Quinoa-BC can reduce the bioaccumulation of Cd in cacao leaves in both moderately acidic and slightly alkaline soils with varying concentrations of the contaminant. On average, leaf Cd was 1.8-and 2.5-fold lower when Quinoa-BC at 2% was applied on slightly alkaline and moderately acidic soils, respectively. For Figure 2. Box plot for reduction factor (Rf) in coffee-, inoculated-, quinoa-, and palm-biochar amended soils at 1 and 2% rate at 130 days after greenhouse experiment. Non-spiked (a) and Cd-spiked (b) slightly alkaline soils. Non-spiked (c) and Cd-spiked (d) moderately acidic soils. Red line represents 1, values greater than 1 have a positive effect, i.e., lower Cd in plants, whereas values lower than 1 is the opposite (Vanderschueren et al., 2021). Values represent the mean (n ¼ 3). Different letters indicate statistical significance at P < 0.05 (Tukey's test). the acidic soil, the reduction of leaf Cd was higher than the one reported in Ramtahal et al. (2019) (2.5 vs 0.9-fold) with the application of commercial biochar (652 kg ha À1 ) after 120 days. Cadmium concentration in leaves were negatively correlated (P < 0.01) with soil properties, i.e., soil pH, CECe, and EC. The increase of soil pH in Quinoa-BC amended soils decreased bioavailable Cd, as discussed above, which consequently reduced Cd uptake in cacao plants. Aside from Quinoa, Only Palm-BC was able to reduce Cd bioaccumulation in leaves, but this was observed only in the Cd-spiked moderately acidic soil. In a batch experiment, with high-soluble Cd concentration, great Cd sorption onto the Palm-BC surface was observed (L opez et al., 2020). Thus, the Palm-BC at a 2% rate, in addition to Quinoa-BC, could potentially alleviate excessive bioavailable Cd in moderately acidic soils with total soil Cd higher than 1.60 mg kg À1 . Cadmium concentration in leaves can be considered as a proxy for Cd bioaccumulation in cacao beans since these two tissues are highly correlated. For instance, correlation coefficients of 0.74 (Ar evalo-Gardini et al., 2017) and 0.91 (Argüello et al., 2019) have been reported. Our results showed that Quinoa-BC effectively reduced leave-Cd in all soil conditions thus this amendment should be tested in field remediation trials.
To understand the effectiveness of biochars in lowering plant-Cd, the Rf values were computed. A Rf of 1.5 could have significant impact on cacao plantations affected by elevated Cd levels (Vanderschueren et al., 2021). With the application of Quinoa-BC an average Rf of 1.6 and 2.0 was observed in slightly alkaline and moderately acidic soil, respectively. In our study, the Rf found in moderately acidic soil is higher than the one found by Ramtahal et al. (2019) (2.0 vs 1.5) with the application of commercial biochar (652 kg ha À1 ) after 120 days. Azar et al. (2019) found a similar Rf (2.1) in wheat and rice plants using rice husk biochar at a 2% dose. A little higher Rf (3.7) was attained in a neutral-pH soil with the application of Miscanthus x giganteus biochar. To achieve such reduction, authors applied a much higher dose (10%) as compared to the one applied in our study. Indeed, at the lowest dose (1%), authors reported that Cd transfer was enhanced (Rf < 1) (Houben et al., 2013). This is in line with our findings, some biochars (e.g. Inoculated-BC) at low dose increased Cd transfer from soil to plants likely due to binding in low-weight organic molecules (O'Connor et al., 2018). Application of quinoa biochar at 2% reduced Cd uptake by plants by approximately 1.8-fold compared to the control. This reduction in the mass of Cd (uptake amount) implies a lower amount of contaminant entering the food chain, due to the mass of Cd immobilized in the biochar-amended soil.
Selecting a realistic rate is critical for studying the effect of biochar for Cd immobilization. It has been widely documented that high rates, e.g. 5% w:w, can reduce Cd transfer from soil to plant parts (Zhang et al., 2016;Dai et al., 2017). However, such application rates are not affordable under agricultural field conditions (Dai et al., 2017). In our study, Quinoa-BC at a realistic dose of to 40 t ha À1 (equivalent to 2%) showed to be a potential amendment to reduce Cd bioaccumulation in cacao plants in field studies. Biochar derived from quinoa has the characteristics to alter soils properties which reduces Cd accumulation in cacao plants. Figure S4 shows the relationship between biochar ash, basic cations content, EC, and Rf of the studied biochars and other biochars reported in the literature (Houben et al., 2013;Mohamed et al., 2018;Xiao et al., 2019). The results indicated that these properties of biochars are significantly correlated (r ¼ 0.66 to 0.85) with their effectiveness in reducing Cd bioaccumulation in plants. Likewise, Person's correlations showed a significant (P < 0.01) and negative relationship between biochar properties such as ash and basic cations content and CaCl 2 -, NH 4 OAc-, and HCl-Cd concentration (r ¼ -0.50 to -0.96) (Table S7). Previous research reported that the ash content of biochars plays an important role in the sorption of metals Yuan et al., 2019). In line with these results, biochars basic cations concentration showed to be a useful predictor of Cd immobilization in aqueous solutions (L opez et al., 2020). The characteristics of the biochars are critical for selecting an amendment suitable for field trials. It is suggested that ash content could be used as a guideline in the identification of potential biochar for Cd immobilization in soils. Figure 3. Relationships between soil pH and CaCl 2 extractable Cd. The data corresponds to the soil pH of coffee-, inoculated-, quinoa-, and palm-biochar amended soils at 1 and 2% rate at 130 days after greenhouse experiment and values for CaCl 2 extractable Cd. Data include values for non-spiked and Cdspiked moderately acidic and Cd-spiked slightly alkaline soils. Figure does not take into account the CaCl 2 -Cd values of the non-spiked slightly alkaline soil. Table 4. Pearson's correlation analysis between soil properties and CaCl 2 -, NH 4 OAc-, and HCl-Cd. The data corresponds to the CaCl 2 -, NH 4 OAc-, and HCl-Cd values and soil properties of non-spiked and Cd-spiked moderately acidic soil and Cd-spiked slightly alkaline soil amended with coffee-, inoculated-, quinoa-, and palm-biochar at 1 and 2% rate at 130 days after greenhouse experiment. Statistical significance is represented by *P < 0.05, **P < 0.01. y NH 4 OAc-exchangeable basic cations.

Biochar application enhances nutrient levels in cacao plants
Aside from potential toxic elements (e.g. Cd) immobilization, biochars could improve soil fertility and plant nutrition (El-Naggar et al., 2019). The uptake of these essential plant nutrients can boost cacao yields. Remediation of Cd contaminated soils can be more sustainable if the amendments can provide essential plant nutrients and impact yields (Hamid et al., 2019). Our results indicated that the application of Quinoa-BC reduced plant-Cd and improved P and K uptake simultaneously (Table S8). Significant increase in plant-P and -K concentration in moderately acidic soil indicates the nutrient supply by the application of Quinoa-BC which could reduce the application rate of P and K fertilizers. Similarly, in slightly alkaline soils, an increase in plant-K was observed (Table S8). Similar increases in plant nutrients have been reported in other studies as a side effect of biochar application. Bashir et al. (2018) found that the addition of biochar at 3% significantly increased N, P and K concentration in Chinese cabbage (Brassica rapa). The application of bamboo chips and rice straw biochar at 2.5% significantly increased P and K content in Brassica chinensis. Remediation of Cd contaminated soils can be more sustainable if the suitable amendments can also provide plant nutrients and potentially increase yields (Hamid et al., 2019).

Conclusions
There is an urgent need to find agronomic countermeasures to alleviate Cd bioavailability in cacao growing areas. The incorporation of biochar derived from quinoa residues at a 2%, increased soil pH and effective cation exchange capacity which are key soil properties governing bioavailable Cd. Quinoa derived biochar display a lower bioavailable soil-Cd due to mechanisms such as co-precipitation and ion exchange decreasing leaf-Cd. Additionally, application of quinoa derived biochar provided plant mineral nutrients (P and K), which could increase productivity to offset mitigation costs. The selection of biochar based on its physicochemical properties, prior to field experiments, is essential to increase the success rate of Cd immobilization treatments on cacao farms. With our results it is possible to conclude that the ash and basic cations content, and functional groups could be good indicators for the selecting or engineer biochars to reduce Cd bioavailability in areas with diverse soil conditions. The results can be applied as a basis for the selection of materials for trials on cacao farms. A field experiment with quinoa derived biochar at a comparative rate, i.e., 40 Mg ha À1 , should be implemented to evaluate the findings.

Author contribution statement
Juli an E. L opez: Conceived and designed the experiments; Performed the experiments, analyzed and interpreted the data; Wrote the paper.
Catalina Arroyave & Adriana Aristiz abal: Conceived and designed the experiments; Analyzed and interpreted the data.
Eduardo Chavez: Conceived and designed the experiments; Analyzed and interpreted the data; Contributed reagents, materials, analysis tools or data.
Catalina Arroyave was supported by Patrimonio Aut onomo Fondo Nacional de Financiamiento para la Ciencia, la Tecnología y la Innovaci on Francisco Jos e de Caldas [120680863411].
Juli an E. L opez was supported by Scholarship [201706101].

Data availability statement
Data included in article/supp. material/referenced in article. . The data corresponds to the CaCl 2 , NH 4 OAc, and HCl extractable Cd values and leaf-Cd concentration for plants of cacao cv CCN-51 planted on non-spiked and Cd-spiked moderately acidic soil and Cd-spiked slightly alkaline soil amended with coffee-, inoculated-, quinoa-, and palm-biochar at 1 and 2% rate at 130 days after greenhouse experiment. Figure does not take into account the CaCl 2 -Cd values of the non-spiked slightly alkaline soil.

Declaration of interest's statement
The authors declare no conflict of interest.

Additional information
Supplementary content related to this article has been published online at https://doi.org/10.1016/j.heliyon.2022.e09790.